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Abstract 

We find sets of solutions to the generalized spheroidal wave equation (GSWE) or, equivalently, 
to the confluent Heun equation. Each set is constituted by three solutions, one given by a 
series of ascending powers of the independent variable, and the others by series of regular and 
irregular confluent hypergeometric functions. For a fixed set, the solutions converge over dif- 
. ferent regions of the complex plane but present series coefficients proportional to each other. 

', These solutions for the GSWE afford solutions to a double-confluent Heun equation by a 

Q ■ taking-limit process due to Leaver. Another procedure, called Whittaker-Ince limit, provides 

2^ . solutions in series of powers and Bessel functions for two other equations with a different type 

of singularity at infinity. In addition, new solutions are obtained for the Whittaker-Hill and 
CN ' Mathieu equations by considering these as special cases of both the confluent and double- 

confluent Heun equations. In particular, we find that each of the Lindmann-Stieltjes solutions 
psj i for the Mathieu equation is associated with two expansions in series of Bessel functions. We 

' also discuss a set of solutions in series of hypergeometric and confluent hypergeometric func- 

■ tions for the GSWE and use their Leaver limits to obtain infinite-series solutions for the 

, Schrodinger equation with an asymmetric double-Morse potential. Finally, the possibility of 

extending the solutions of the GSWE to the general Heun equation is briefly discussed. 

o 

^ ; I. Introduction 



We study sets of solutions for a generalized spheroidal wave equation (GSWE) aiming at 
obtaining solutions for other differential equations which are related to the GSWE by taking-limit 
processes. Although the three solutions which constitute a fixed set converge over different regions 
of the complex plane, they have series coefficients proportional to one another what means that, 
- under certain circumstances, the three solutions imply a unique characteristic equation resulting 
^> ! from the three-term recurrence relations for the coefficients. Solutions for a double-confluent Heun 
\ equation (DCHE) are generated from the ones for the GSWE by using of a limit process devised by 
Leaver [1] and, therefore, without solving directly the DCHE. On the other hand, by applying to 
the GSWE and DCHE another limit due to Whittaker and Ince [2l[3lll], we obtain two differential 
equations with a different type of singularity at infinity whose solutions are obtained from the ones 
for GSWE and DCHE by taking the Whittaker-Ince limit, again without integrating directly the 
equations. Thus, by means of the Leaver and Whittaker-Ince limits, the GSWE afford solutions 
for three other different equations. 

Furthermore, new solutions for the Whittaker-Hill and Mathieu equations are obtained from 
the fact that these equations are particular cases of both the GSWE and the DCHE [5]. From 
one side, this gives two types of solutions for the Whittaker-Hill equation (WHE) and, from the 
other side, three types for the Mathieu equation in virtue that this one is also an instance of the 
Whittaker-Ince limit of the GSWE. 

In summary, solutions for the GSWE generate solutions for five other equations which in most 
of the cases have been studied separately and without any connection with the GSWE, namely: 
DCHE, Ince's limits of the GSWE and DCHE, and the Mathieu and Whittaker-Hill equations. 
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Some solutions we shall find are already known, but here they are associated with two other 
solutions possessing different mathematical properties, since we are dealing with sets containing 
three solutions, as mentioned in the first paragraph (the known solutions result from a known 
solution for the GSWE given by a series of ascending powers of the independent variable). 

In addition to the solutions for the GSWE and its limiting and particular cases, we consider 
two physical problems. One is given by the Schrodinger equation for asymmetric double-Morse 
potentials concerning quantum spin systems [6l [7] . For quasi-exactly solvable potentials [H [9] , an 
earlier attempt to solve this equation had been unsuccessful [10]. Here we show that the solution 
of this problem requires the use of two different types of solutions for the DCHE, both of them 
discussed in this article. The other problem is given by the Teukolsky equations for gravitational 
backgrounds of black holes, which afford a connection among three equations: general Heun 
equation, GSWE and DCDE [H], [T2]. However, this problem is regarded only at the extent it 
helps to find solutions for the Heun equation and its confluent cases. 

After these general considerations, let us examine each of these six equations and how they are 
connected, discuss some features of the solutions and outline the structure of the paper. 

The generalized spheroidal wave equation (GSWE). For the GSWE we adopt the form used by 
Leaver [T], 

z{z -^o)^ + (5i + B2z)— + [53 - 2r]u;{z - zq) + uo\{z - zq)] f/ = 0, (c^ ^ 0) (1) 

where Bi, rj and uj are constants and z = and z = zq are regular singular points with indicial 
exponents (0, 1 + Bi/zq) and (0, 1 — B2 — Bi/zq), respectively, that is. 



lim U{z) ~ 1 or lim U{z) z , 



2-^0 2-»0 

vl-B2-(Bi/2o) 



l+(Bi/2o) 

(2) 



lim U{z) ~ 1 or lim U (z) ^ {z — zq 

z^zo z—*zo 

At 2; = 00, which is an irregular singularity, the behavior inferred from the normal Thome solutions 
[I1[I3] is 

lim U{z) ~ e^'^'z^'''-^^'/^\ (3) 
2^00 

Eq. ([1]) is in a form appropriate for the study of the Teukolsky equations for the Kerr metric 
since, in the upper limit for the rotation parameter, zq vanishes and the equation reduces to an 
DCHE. In fact, Eq. ([1]) is equivalent to the confluent Heun equation [5l[Tll[T5] and, when rj = 0, 
it reduces to the so-called ordinary spheroidal wave equation [T6] . 

The double-confluent Heun equation (DCHE). By setting 2:0 = in Eq. ([T]), Leaver obtained a 
DCHE with five parameters fl], namely. 



z 



+ (5i + B2Z) -^+{B3- 2VUJZ + u'z'-) ?7 = 0, (Si ^ 0, ^ 0) , (4) 



where z = and ^ = 00 are both irregular singularities. At 2; = 00 the behavior of the solutions 
is again given by Eq. ([3]), while at 2; = the normal Thome solutions affords 

limf/(z) ~ 1, or \imU{z) ~ e^'^'z^~^\ (5) 

2— >0 2—^0 

The values Bi = and/or uj = (for t] fixed) are excluded because: if i?i = = we have the 
Euler equation; if i?i = and 7^ or if i?i 7^ and lj = 0, the equation degenerates into a 
confluent hypergeometric equation (Appendix A of j^). 
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Here the Leaver limit 2:0 ^ is used to obtain solutions for the DCHE (jlj) from solutions of the 
GSWE. Furthermore, the Leaver form (jl]) for the DCHE has also the advantage of admitting the 
Whittaker-Ince limit, in opposition to other forms which have only four parameters, as in Refs. 

13 [a [n]. 

The Whittaker-Ince limit of the GSWE. We define the Whittaker-Ince limit of the GSWE and 
DCHE as 

— 0, — > 00, such that 2rjuj = —q, (Whittaker-Ince limit) (6) 

where g is a constant. In a previous article we have referred to the above as 'Ince limit', but 
according to Humbert [2], prior to Ince [3], Whittaker had used such procedure to get the Mathieu 
equation from the Whittaker-Hill one. In any case, the limit ([6]) is a generalization of the limit 
used by Whittaker and Ince. 

The Whittaker-Ince limit of the GSWE is 

(fU dU 
z{z - ^0)^ + {Bx + B^z)^ + [53 + q{z - z^)\ t/ = 0, (g ^ 0) (7) 

(if g = this equation can be transformed into a hypergeometric equation) . At the regular singular 
points z = and z = zq the behaviors of the solutions are formally the same as those for the 
GSWE, but at the irregular singular point z = 00, from the subnormal Thome solutions [13] we 
get 

hm U{z) ~ e±2.v/^^(i/4)-(B./2)_ /g) 

2— +00 

in contrast with the behavior ([3]) of original GSWE (normal Thome solutions). 

Mignemi [18] and Malmendier [19] have found a particular case {B2 = 2Bi, zq = —1) of Eq. 
(I7j) for a wave equation resulting from the Laplace-Beltrami operator for a scalar field on the 
Eguchi-Hanson space [20], but they did not succeed in deriving solutions from solutions of the 
GSWE. Here and in Ref. [1] this goal is attained by starting with solutions for the GSWE which 
admit both the Leaver and the Whittaker-Ince limits. One can check that this is not possible for 
some other solutions as the Hylleraas and Jaffe ones [1], for example. 

The Whittaker-Ince limit of the DCHE. The Whittaker-Ince limit of the DCHE is the equation 

(fU dU 
z^-^ + (5i + B2Z)— + (53 + qz)U = 0, (g ^ 0, 5i ^ 0) (9) 

where 2; = and 2; = 00 are irregular singularities as in the DCHE (jl]). At the irregular singular 
point 2; = the solutions behave again according to Eq. dSD, but at infinity they admit the 
behavior given by the subnormal Thome solutions, Eq. ([8|). 

Eq. ([9|) is not obtainable from a DCHE with four parameters. In fact, the scale transformation 
in the variable z used to pass from five into four parameters [5l HI] becomes meaningless when 
= in Eq. (jlj). So, the Leaver form (jlj) is valuable since it leads to solutions for Eq. (j9j) 
which rules, for example, the radial dependence of the Schrodinger equation for the scattering of 
low-energy ions by polarizable neutral atoms with an induced quadrupole momentum [1]. 

In Eq. (jlj) g = and/or Si = were excluded because in these cases the equation can be 
transformed into a confluent hypergeometric equation if g = and -Bi 7^ 0, into a Bessel equation 
if g 7^ and -Bi = 0, and into an Euler equation if g = _Bi = [1] . 

The previous equations - (jlj), (jlj), (jZj) and (jHj) - can also be distinguished by the rank or species 
of the singularity at z = cxd [5], jH]. For the GSWE (jlj) and the DCHE (jlj) the rank is 1 {species 
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2), while for the Whittaker-Ince hmits ([7]) and the rank is 1/2 {species 1). However, here this 
classification is not relevant. To get the solutions, more important are the concepts of Leaver and 
Whittaker-Ince limits, and the information on whether the solutions exhibit the behavior predicted 
by either the normal or the subnormal Thome solutions. 

Notice that Eq. ([9]) is also the Leaver limit of Eq. (I7j) and, consequently, we have two ways to 
generate its solutions, as indicated in the following scheme: 



GSWE (1] 



Leaver limit 



DCHE (4) 



i 

Whittaker-Ince limit 

i 



i 

Whittaker-Ince limit 

i 



(10) 



Eq. (7) 



Leaver limit 



Eq. (9) 



A preliminary study of solutions for the Heun equation indicates that above diagram can be 
improved by adding the connection Heun equation — > GSWE (see Appendix B). 

Whittaker-Hill and Mathieu equations. The Whittaker-Hill and Mathieu equations are partic- 
ular cases of both the GSWE ([1]) and the DCHE (jlj), but the Mathieu equation is also a special 
instance of the Whittaker-Ince limit of the GSWE. Thence, there are two types of solutions for 
the Whittaker-Hill equation (WHE) and three types for the Mathieu equation, as stated above. 
This result is due to Decarreau, Maroni and Robert [5] but it is necessary to study the properties 
of the solutions obtained by each of these approaches, by writing such solutions explicitly and 
comparing them. 

To avoid confusion in the notation, we write the Whittaker-Hill and Mathieu equations in a 
form independent of the parameters Bi, r) or u which appear in the other equations. We choose 
the trigonometric (hyperbolic) form because this is useful to investigate the parity and periodicity 
properties of such solutions. For the WHE we employ the form [Sj [2T] 



du^ 



?9 - -^^ - (p + 1)^cos(2ku) + cos(4ku) 



iy = 0, (WHE) 



(11) 



where we have written i?, instead of the usual i], because rj has already appeared in the GSWE 
and DCHE. If m is a real variable, this equation represents the WHE when k = 1 and the modified 
WHE when k = i. The WHE ^ is converted into a GSWE or DCHE in Sees. II B or HI B, 
respectively. 

For the Mathieu equation we use the form [22] 



(PW_ 
du^ 



+ cr^[a - 2/c^ cos{2au)]W = 0, (Mathieu equation) 



(12) 



where a = 1 or a = i for the Mathieu or modified Mathieu equation, respectively. Note that there 
are two ways to get Eq. f|T2l) from the WHE f|TT]) . First, for p = — 1 and k = a/2, the WHE 
(fTTl) reduces to the Mathieu equation ( fT2l) . Second, the Mathieu equation (fT2l) also follows from 
the WHE ( ITTi) by means of the original Whittaker-Ince limit [21 [3l |23]: ^ — >^ 0, p — oo so that 



= 2k' 



K 



a and ^ = a. This second way is equivalent to consider the Mathieu equation 



as a particular case of the Whittaker-Ince limit ([7]) of the GSWE. The Mathieu equation (fT2|) is 
transformed into a GSWE, a DCHE or the Ince limit of the GSWE in Sees. II C, HI C or IV B, 
respectively. 

In this paper, starting from sets formed by three solutions for the GSWE, we follow the scheme 
(llOp to generate solutions for the other equations. In each set, one solution is given by a Barber- 
Hasse expansion in power series [U |2l] while the others are given by expansions in series of regular 
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and irregular hypergeometric functions. In a fixed set, the solutions converge over different regions 
but are collected together because their series coefficients are proportional to each other. The 
construction of sets of solutions with these features is in part suggested by study of the Mathieu 
equations [22j, but it follows as well from the study of wave equations in general relativity where 
sometimes is necessary to match solutions convergent over different domains [U Uni EHl EHl [271 EH] ■ 

In the context of the GSWE, a initial set, as the one given in Eq. f l33al) . leads to the others 
by means of transformation rules coming from variable substitutions which preserve the form of 
the GSWE but modify their parameters and/or arguments. The transformation rules for the 
GSWE, written in Eqs. (I13ap . allow us to generate solutions with all the expected behaviors at 
the singular points. These features are transferred to the solutions of the limiting equations and 
their particular cases. This is illustrated by the results of Sections II B and IV B, where we find 
that each set of solutions for the Whittaker-Hill and Mathieu equations, respectively, possesses 
different properties of parity and /or periodicity. 

On the other hand, the convergence of the series depends on whether the parameters of the 
differential equation do or do not satisfy a relationship, known as characteristic equation, which 
results from the three-term recurrence relations for the series coefficients. In fact, there are two 
cases to be considered. First, if there is some arbitrary constant in the differential equation, the 
characteristic equation can be satisfied by using it to determine that constant in terms of the 
others and in such case the series converges. If there is no disposable constant, the characteristic 
equation in general cannot be satisfied and then the series diverges. In this second case, the con- 
vergence is assured by inserting into certain solutions a parameter z/, called phase or characteristic 
parameter, which plays the role of the arbitrary constant and must be adjusted so as to validate 
the characteristic equation. Excepting section VI, all the other sections of the present article deal 
with solutions without characteristic parameter and, consequently, assume the existence of an 
arbitrary constant in the differential equations. 

In section II A, we study the Barber-Hasse solutions and the expansions in series of confluent 
hypergeometric functions for the GSWE. In section II B these solutions are particularized for the 
Whittaker-Hill equation and in Sec. II C for the Mathieu equation. In each set of solutions for the 
WHE, the solution resulting from the Barber-Hasse one turns out to be a solution already found 
by Arscott [21j, but the others seem to be new. 

In Sec. HI A the Leaver limit is used to derive solutions for the DCHE from solutions for 
the GSWE. The solutions ensuing from the Barber-Hasse solutions include as particular cases the 
ones which have been investigated by Decareau, Maroni and Robert [Hj and also by Schmidt and 
Wolf [T7j for a DCHE with four parameters, but to each of these now correspond two expansions 
in series of confluent hypergeometric functions. These solutions are used to find finite-series 
solutions for the Schrodinger equation with an asymmetric double-Morse potential. There are 
infinite-series solutions but these are not regular at the singular points of the equation. The WHE 
and the Mathieu equations as particular cases of the DCHE are discussed in Sees. HI B and HI 
C, respectively. 

In Sec. IV we deal with the Whittaker-Ince limit of the GSWE and find that the expansions 
in series of confluent hypergeometric functions reduce to expansions in series of Bessel functions. 
In the special case of the Mathieu equation, we recover some solutions found by Lindmann and 
Stieltjes ED], each one coupled with two expansions in series of Bessel functions. Similarly, in 
Sec. V we establish the Ince limits for the solutions of the DCHE. 

In Sec. VI we consider some solutions with a phase parameter for the case in which there is 
no free constant in the GSWE. Each set of solutions for the GSWE contains two expansions in 
series of hypergeometric functions and two in series of confluent hypergeometric functions. These 
solutions also admit of both the Leaver and Whittaker-Ince limits. The solutions for the DCHE 
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are used to find regular infinite-series solutions for the Schrodinger equation with the asymmetric 
double-Morse potential mentioned above. For this it is necessary to match two solutions converging 
over different ranges of the independent variable but having a common region of convergence. The 
results is that, for this quasi-exactly potential, the entire spectrum of energies may be determined 
by solving a characteristic equation given by the sum of two infinite continued fractions. 

Section VII contains some concluding remarks, while Appendix A presents an alternative 
derivation of the expansions in series of Bessel functions for the Whittaker-Ince limit of the GSWE 
and Appendix B discusses some connections between solutions of the general Heun equation and 
generahzed spheroidal equation. 



II. Solutions to the generalized spheroidal wave equation 

In this section, first we write down transformation rules for the GSWE and recall some basic 
properties of the three-term recurrence relations for the series coefficients (not restricted to the 
GSWE). After this we set the solutions for the general case and discuss both the Whittaker-Hill 
and Mathieu equations from the viewpoint of the GSWE. 

In the first place, if U{z) = U{Bi, B2, B^, zo,uj,ri] z) denotes one solution of the GSWE, the 
transformation rules Ti, T2, T3 and T4 are given by [H [3Tj 

T,U{z) = z'+^-/^'^U{C,,C2,Cs;zo,u,v;z), zo^O, 

T^Uiz) = {z- ^o)^-^^-^^/"«f/(Bi, D;; zo, to, V, ^07^ 0, 

(13a) 

TsUiz) = U{Bi, B2, B3; Zq, -uj, -T]; z), Vzq, 

T^U^z) = U{-Bi - B2Zo,B2,B3, + 2r]ujZo; zq, -uj,t]; zq - z), Wzq 

where 

C, = -B,-2zo, C2 = 2 + B2 + '-^, C3 = 53+(l + f)(52 + f), 

(13b) 

D2 = 2-B2-'-^, ^3 = 53 + f (f +52-1). 

These rules can be verified by substitutions of variables. In T3 it is assumed that we must change 
the sign of {f],uj) only where these quantities appear explicitly in the solution U{z), keeping the 
expressions for the other parameters unchanged even if they depend on rj and u as in the Teukolsky 
equations of the relativistic astrophysics [I],[32]. For brevity, we use only Ti and T2. The rule T4 
interchanges the position of the regular singular points {z = zq ^ z = 0) and may be used to get 
solutions with the appropriate behavior at z = 0. 

On the other hand, if 6„ denotes the series coefficients of a solution, then the three-term 
recurrence relations (for solutions without phase parameter) have the form 

"0^1 + Pobo = 0, anbn+l + PnK + InK-l = (^ > 1), (14) 

where a^, /3„ and 7,^ are constants depending on the parameters of the differential equation. 
These recurrence relations constitute a infinite system of homogeneous linear equations which can 
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be written as 



( I3q aQ ■ ■ 

7i (3i "1 

72 /?2 







a2 



'~in Pn Clji 



v 



hi 
b2 



(15) 



Then, to have nontrivial solutions for the (Hill) determinant of the above infinite tridiagonal 
matrix must vanish. This is possible only if there is some arbitrary parameter in the differential 
equation, that is, in the elements of the matrix. Equivalently, the recurrence relations lead to 
characteristic equation in terms of the infinite continued fraction [1] 



l3o 



ao7i '^172 ^273 
A- /32- /^s- 



(16) 



So, instead of a vanishing determinant, we must satisfy the characteristic equation and this also 
assures the series convergence by means of some version of a Perron theorem |33j . 

Under certain conditions, three-term recurrence relations like (1141) give finite-series solutions 
which are also called quasi-polynomial solutions, Heun polynomials or quasi-algebraic solutions. 
The condition sufficient for having terminating series is 7„ = for some n = N =positive integer 
[M] . in which case the series presents terms, that is, < n < iV — 1 . In effect, if 77V = we can 
choose the parameters of the equation so that = 0: then, the relations (|T^ imply that bn = 
for any n > N. In such event, in general there are infinite series as well since the transformation 
rules generate solutions with other recurrence relations. In addition, if a„ = for some n = N, 
the series begins at n = + 1, but in this case one may take n = m + N + 1 and relabel the series 
coefficients in order to obtain a series beginning at m = 0. 

There is another remark which will be useful in Sees. II A and IV A. Suppose a second solution 
with coefficients c„ having the recurrence relations 



ttoCl + PoCq = 0, a„Cn+l + PnCn + lnCn-1 = (?2 > 1), 

where /3„ is the same as in Eq. iHM . Then, if 



(17a) 



(17b) 



it is obvious from Eq. (ITB]) that, for infinite series, both solutions have the same characteristic 
equation provided that the sums begin at n = in both series. 



A. GSWE and the Barber-Hasse solutions 

Now we construct the initial set containing the solutions ([/{*, f/f°, Ui). Ui is the Barber-Hasse 
expansion in series of ascending power oi z — zq which converges for any finite z and was originally 
proposed to solve an angular GSWE in which < 2; < 2:9 tlj- ^1° is given by an expansion in series 
of irregular confluent hypergeometric functions \l/(a, 6; y) (or confluent hypergeometric function of 
second kind), converges for \z\ > \zq\ and, thus, is suitable to solve a radial GSWE {z > zq) except 
at z = Zq. Ui is an expansion in series of regular confluent hypergeometric functions $(a, b; y) (or 
confluent hypergeometric function of first kind) which converges for any value of z but, in general. 
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its behavior at infinity is not the one predicted by the Thome solutions. From the first set of 
solutions we obtain three other sets by using the transformation rules Ti and T2. 

To obtain the Barber-Hasse solution f/° we write U{z) = Ui{z) = e^^''^ f{z). Then, the GSWE 
([1]) leads to 

(Pf df 
z{z -zo)-^ + [Bi + B2Z + 2iujz{z - zo)] ^ + [^3 + iujBi + iujB2Z - 2ujr]{z - zq)\ / = 0. (18) 

Now, by expanding f{z) in a Frobenius series corresponding to the exponent about z = zq, i. e.. 



Z - Zo 



U',iz) = e-J2b^:\z- 



(19) 



n=0 



n=0 



we find that bn^ satisfy the recurrence relations {b]^{ = 0) 



(1) 



n{n + B2-l + 2iujzo) + iuzo [B2 + ^] + B3 



20 



+2iuj {n-l + i7] + f) b^^\ = 0. 



(20) 



This solution can also be obtained by applying the rule T4 to the Barber-Hasse expansion given in 
Ref. [Jj. It converges for any finite z provided that the characteristic equation is fulfilled [1] |2^. 

To find the solutions in series of confluent hypergeometric functions, we carry out the substi- 
tutions 



t/(^) = e-^F(y), y = -2zu;z, 
which transform the GSWE ([T]) into 



(21a) 



{y + 2iLJZo) 



dF 



- y- 



dy'^ dy 



^ dF ^. „ dF 

+ B2y—-2iujBi— + 
dy dy 



i?3 + iujBi + 2rnjjZQ — [ir] + 



Bo 



y 



F = 0. 



(21b) 



Then, first we expand F{y) in series of the irregular functions \l/(ci, b; y) [35] according to 

n=0 ^ ^ ri=0 

f B \ 

V^{z) = e'"^' ci^)^ n + 27/ + n + B2; -2iujz , 

where '^niy) = ^ {n + irj + B2/2,n + B2]y). These '^niy) satisfy 



(22) 



-(n + B2 



,d'fu{y) 



dy 
y dy 



2J dy 

(n + iri + f) ^Hn+iiy), 



+ {n + tr] + f) ^„(y). 



(23) 



{l-n-B2 + y)^n{y)-'^n-i{y), 
{n + zr] + f) y^n+M + (1 - n - ^2 + y)^„(i/) - ^n-i(z/) = 0, 
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where the first relation is the confluent hypergeometric equation for \l/„(?/) and the others follow 
from the properties of the irregular confluent hypergeometric functions [35] . Thus, inserting the 
previous expression for F{y) into Eq. fl21bl) and using Eqs. fl23l) . we obtain 



n=l n=0 



2:0 



n + iri + ^] ^„+i(y) 



0,(24) 



where /3n is the coefficient of bn in the recurrence relation for the solution U^, namely. 



/?(^) = n{n + B2-l + 2iuzo) + iuzo ( ^2 + — ) + ^3. 



Putting m = n — 1 and m = n + 1 in the first and last terms of Eq. (1241) . respectively, we find 



E 

m=l 



(m + l)c^?,, + 4^)c« + 2^a;2;o i m - 1 + + — ] i m - 1 + tr] + ^ ) c^^l. 



Bi 



Bn 



2 \ ^(1) 



= 0. 



By equating to zero the coefficients of each independent '^n{y) we find the three-term recurrence 
relations {c^^\ = 0) 



n + l)ci'ii + /3(^)cW + 2iuJZo { n - 1 + B^ + — ] { n - 1 + tr] + ) c^^'i, = 0, n > 0. (25) 



Bn 



2 \ ^(1) 



Zo 

From the behavior of \I^(a, 6; y), when h ^ 00 while 6 — a and ?/ remain bounded [35], we get 

lim H?±iM = J_. 

n^oo \l/„(?/) 2iujz 

Besides this, using a Perron-Kreuser theorem [33,J we find that the minimal solution of the recur- 
rence relations (l25l) satisfies 



lim 



.(1) 



Then, combining these two limits, we have 

lim 



-2iujZo, if Zo 7^ 0; 

-2iujBi/n, if zq = 0. 

-2:0/2:, if 7^ 0; 

(n^;), if Zq = 0. 



(26) 



Therefore, by the ratio test, the series in U^{z) converges for any \z\ > \zq\ on the condition that 
the characteristic equation is satisfied. The behavior of U°°{z) when z —>■ 00 results from the 
relation |35l 



lim ^(a, b; y) ~ + 0{\y\-% -— < argy < - — , 



y-*oo 



(27) 



which implies 



lim U^{z) ~ e"^'z 



iujz -iri-{B2/2) 



3tt Stt 
— < axg{-2iujz) < 



as in a normal Thome solution. 

Now we expand F{y) in series of regular confluent hypergeometric functions ^{a,b;y). If we 
define the function as 

^n{y) = ^In + B^f (^^ + + Y' + ^ ) ' ^ = "^'"""^ ^^^^ 

and use some difference and differential properties of $(a, 6; y) [36], we find that ^n{z) also satisfy 
the relations (123|) . Therefore, substituting for '^niy) in Eq- (|22|) we get the solution 



n=0 n=0 ^ ^ 



(29) 



where the coefficients ch^ are the same which appear in the solution U^{z). In order to show that 
the series in Ui{z) converges for any z, we use the behavior of $(a, b; y) when b —>■ oo while b — a 
and y remain bounded [35j, namely. 



lim $(a, 6; y) = e^^fl + 0{\b\-% {b - a = finite) 

6— >oo 

which leads to 

lim *?±iM . 

Thus, using again the minimal solution of the recurrence relations for Cn we get 



Therefore, the ratio test implies that the series in the solution Ui{z) converges for any value of z. 
However, the behavior of Ui{z) when 2; ^ oo in general does not coincide with the one inferred 
from the Thome solutions, since we have [36] 

r ig)ev-ni+ 0(12/1-1)], (%>o) 

lim $(a, b;y)= { (30) 

i fiS)(-?/)-1i + o(li/|-^)], (%<o). 

Thus, only if 3??/ < the asymptotic behavior of Ui has the same form as . If 3ft?/ = 0, the 
limit of ^(a, 6; y) is a combination of the two expressions given on the right-hand side of the 
previous expression [35]. On the other hand, although the function $(a,6;y) is not defined if 
6 = 0, —1, —2, ■ ■ ■ , the function <l>(a, 6; y) is defined since [35] 

lini ^i^lM = M!!iy-$(a + m,l + m;y), m = l,2,3,---, (31) 

(o)m = fl(fl + 1) ■ ■ ■ (a + m — 1) for m > 1. Therefore, the expansion Ui{z) is defined even if B2 is 
zero or a negative integer. 

In the following we collect the three solutions in the first set of solutions, and this gives three 
other sets by the use of the transformation rules Ti and T2 according to 

{ui f/r , f/i) ^ {ui f/r , U2) ^ {ui u^, f/3) ^ {ui ur, u,) ^ {ui u^, u,) . (32) 
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First set. This set admits both the Leaver and Ince hmits. If B2 + Bi/zq = 1, it is equal to the 
fourth set. 



00 



n=0 
00 



Setting b^^l = and cl[ = 0, the three-term recurrence relations for the coefficients are given by 



/ B \ 

U^{z) = e^^^^ci^)^ U + ir] + ^,n + B2;-2iuz] , 

n=0 ^ ^ 

Ui{z) = e''^' in + ir] + ^,n + B2] -li^z J . 

n=0 ^ ' 

.(1) 



(33a) 



( n + ^2 + ^ ) + \)h^ll^ + + ^2 - 1 + 2%UJZ^) + %UJZ^ ( ^2 + ) + 53 

+2itu - 1 + zr/ + ^) = 0, 



(33b) 



and 



(ra + 52 - 1 + 2zco'Zo) + ^^^0 ( -B2 + -j^ ) + -B3 



22Co'2;o ( ^ — 1 + -B2 



5i 

2^0 



n — 1 + i?7 + 



5, 



2 '\ .(1) 



0. 



The coefficients 6n ^ and Cn '' are connected by 



(1) 



.(1) 



rfn + S2 + — ) 6(1), if i?2 + — 7^0,-1,-2, 

V ^0 / 2^0 



(33c) 



(34) 



apart from a multiplicative constant independent of n. The restrictions on B2 + {Bi/zq) assure 
that the sum in the three solutions begins at n = and, thus, they have the same characteristic 
equation for the reason given in the paragraph containing Eqs. fll7ap and (]17bp . In addition, we 
have finite-series solutions if irj + {B2/2) is zero or negative integer. 

Second set. These solutions admit the Ince limit but not the Leaver limit. If B2 + Bi/zq = 1, 
the second and the third sets are equal to one another. 



f/0(z)=e-V+t5^6f(z-zo)", 

n=0 

^0 ci'^^ (n + l + ir] + — + ^,n + 2 + B2 + 



U^[z) = e'^^^z 



2Bi 



-2%ujz 



n=0 

00 



JJJz) = e'^^z 



^+^0 V" c^^^^ f n + 1 + ir/ + — + n + 2 + ^2 + 



^0 

25i 



(35a) 



-2iujz 



) 



where the recurrence relations are 



z^ (n + 52 + f ) (n + l)6i'|i + 



n ( n + 1 2iijjZQ -h ^2 + ^ ) + iuz^ ( 5 



5, 



Bi 



5, 



(2) 



+2z^ (^n + zr/ + f + f) 6i'2i=0, 



(35b) 
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and 



n 



(n + 1 + 2^u^Zo + + ^) + ^u^zo (b, + f ) + (l + ^) (b, + f ) + i^a] ci^^ 



+2ia; (^n - 1 + S2 + 
The coefficients are connected by 



n + + ^ + 



.(2) 



Third set. These solutions admit both the Leaver and Ince hmits. 



n=0 

U^{z) ^ e^'^^z^^^ {z - ZQf~^^~^ ^cf-^ in + 2 + irj-^,n + A-B2;-2iu;zj , 



.El 



n=0 

00 



e'-'^'^z' ' ^0 (2; — 2:0) 

where the recurrence relations are given by 

zo(n + 2-B2-^) {n + l)b[ 



Y c^2)$ n + 2 + ir;--^,n + 4-52; -2iujz , 

n=0 ^ ^ / 



(3) 
'n+1 



+ 



n (n + 3 + 2ia;;2o - ^2) + iujzo (2 - B2 - + 2 - B2 + B3 
+2iu (n + l + iri-^) b^^l^ = 



^(3) 
<Jn 



and 



(n+l)c2i + 



n (n + 3 + 2ia;^o - ^2) + icu^o (^2 - - ^ j + 2 - ^2 + 5; 
2ia; + 1 - S2 - ^) (n + 1 + i77 - f ) c^^li = 0. 



In this case, we find 



cf = (zorr(n + 2-S2 



B, 



"jfef , if 52 + ^^2, 3, 4, 



Fourth set. This set admits the Ince limit but not the Leaver limit. 

00 



n=0 

00 



C/f (z) = e'^^z - zo)'-'''-^ V c(f)* L + 1 + i77 - ^ - ^, n + 2 - B 

fr;^ V ^0 2 



U^^z) = e'^^'iz - zo) 



1-B2- 



^0 



n=0 

00 



2 ^; — 2ia;2; 

^0 



n=0 



-Bi -B2 7-) 2Bi \ 

n + 1 + zTj , n + 2 — i32 ; —2zujz j , 

Zo 2 zq J 
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with the recurrence relations 



n(n + l + 2..;^o-52-^)+^^^o(2-i?2-f)+f (i?2 + f -1 



5, 



(39b) 



and 



{n + l)citi + 



n n 



+ l + 2^uz,-B,-'-^)+^uJz,[2-B,-^)+^[B, + ^-l)+B, 
2tuzo (n + 1 - i?2 - f ) (n + ^r^ - f - f ) c^^U = 0- 
This time the connecting formula is 



.(4) 



Bo 



Bi 



if + — 7^2,3,4, 

/ Zo 



(39c) 



(40) 



As the rule T4 changes the position of the regular singular points, it may be used to get solutions 



with the suitable behavior at 2; = 0. For these new sets of solutions, the [/, 



5, 6, 7, 8) again 



converge for any finite value of z, the U°° converge for I2; — 2;o| > \zq\, whereas the Ui converge for 
all values of z. 

The previous convergence properties are transferred to the limiting and particular cases of the 
GSWE. The utility of the convergence of the solutions Ui over the entire complex plane requires 
further investigation. In fact, other solutions converging for any z (also giving by series of regular 
confluent hypergeometric functions) have already been found in Sec. VII of Ref. [1]. In any case, 
in physical problems we have also to consider boundary and regularity conditions, in addition to 
convergence. 

B. Whittaker-Hill equation as a GSWE and Arscott's solutions 

By inserting z = cos^(ku) and W{u) = U{z) into the WHE ffTTl) . we obtain the GSWE 



z{z 



+ z 



dU 

dz 



+ 



4^ ^ 4 2 ^ 



z{z 



U = 0. (41) 



From this and Eq. ([T]), we find that the solutions W{u) for the WHE ( ITTi) may be written in terms 
of the solutions U{z) of the GSWE as follows 



W{u) = U{z), z = cos'^{ku) {K = l,i); ^ WHE as 

Zo = 1, fii = -i, B2 = 1, 53 = \[{p+l)^ - zu = l, ir] = Ef J GSWE. 



(42) 



To obtain the solutions for the WHE in terms of z it is sufficient to replace the parameters given 
above into the solutions U{z) of the GSWE, but to get the periodicity and parity properties of 
the solutions we also use the transformation z = cos^(km). 
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First set. Even solutions, period tt ii k — 1 and u —real {ii k — 1 does not converge). 

-1)' 



W^{u) = el -^'('^") J] — A^-— cW sin2"(«xi) , | cos(«ii) | < oo, 

oo 

W^{u) = ei^°^'('^'') J^c^^)* + 1 + |,n + l;-Ccos^(fi;M)) , |cos(fi;M)| > 1, 

n=0 

oo 



(43a) 



n=0 



where the recurrence relations for Cn^ are 



{n + 1)0^1 + 



n{n + + 



2i + pi-^ 



-(1) 



+ (r^ + f)cll = 0. (43b) 



Second set. Even solutions with period 27r if = 1. 

oo 

W^{u) = cos{ku) e^^''^'^^'') ^ 



r[n + (l/2)] 



cos(«;m)| < oo, 



n=0 
oo 



Vl^2~(ii) =cos(k'u) el^°«'('^«)^c(f)* (n+ - + -,n + 2;-^cos2(K'u) ) , |cos(k'u)| > 1, (44a) 

n=o V 2 2 J 

W2{u) = co^{ku) e5^°^'(«")^c(f)$(n+- + |,n + 2;-Ccos2(fi;M) j , Vm 



with the recurrence relations 



(n + l)ci^) + 



n (n + 1 + + 



1 + 2^ + p^-^ 
4 



n 



2 2 



= 0. (44b) 



Third set. Odd solutions, period tt if k = 1. 
Wl{u) = sin(2«:i.) e^-'(-) ^ \ sin'^^iKu), 

n=0 ^ l-^/^^J 



cos(«;ii)| < oo. 



Wi°{u) = sm{2Ku) el^°«'('^") ^c^^)^ (^n + 2 + |, n + 3; -^cos^(fi;M)) , |cos(fi;M)| > 1, (45a) 

n=0 

oo 

Wsiu) = sin(2K'u) ef ^ c^f^l* (^n + 2 + |, n + 3; cos2(«;'u)^ , V-u 



Tt=0 



having the recurrence relations 



(n + l)c£ii + 



n(n + 2 + + H-C + 



+ e + (n + 1 + I) c£l, = 0. (45b) 
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Fourth set. Odd solutions with period 2-7? if k = 1. 



W^{u) = sm{Ku) el^"'^'^-") V / / ^ ^ sm^''{Ku), \ cos(km)| < oo, 

t^o r[n + (3/2)] 

W4°°(m) = sin(ftu) e^=°'^'('^")^c(f)^ ( n + - + |, n + 2; cos2(fi:M) j , |cos(km)| > 1, (46a) 

W4(m) = sm{Ku) e5™«'{''") ^c(f)<l ( n + - + -,n + 2; -^cos2(/tu) J , Vu 



n=0 

where the recurrence relations are 



in + l)c% + 



1 _L r)£ _ ^ 



c':'+^ln+l) (- + ^ + 1)^1 = 0. (46b) 



The solutions corresponding to the opposite choice for the signs of u and rj in Eqs. fH2|) are 
obtained by the substitutions 

e ^ -e, P ^ -P - 2, (47a) 
which take the place of the rule T3. The operation equivalent to the rule T4 is 

u^u + ^, e^-e (47b) 
Similarly, the composition T3T4 corresponds to the symmetry 

71 

u^u + —, p^-p-2, (47c) 
Zk 

of the WHE. In fact, if we apply this last operation to the preceding sets of solutions, we find 
the solutions obtained by Arscott [21], each one accompanied by two expansions in series of 
hypergeometric functions. 

Notice that the use of the transformations Ti and T2 have been essential to generate a subgroup 
containing four sets of solutions which are even or odd and have period tt or 27r. The solutions for 
WHE obtained from the one for DCHE do not possess these properties of parity and periodicity 
(see Sec. Ill B). 



C. Mathieu equation as a GSWE 

Although the Whittaker-Ince limit of the preceding solutions for Whittaker-Hill equation gives 
solutions for the Mathieu equation, such solutions are also particular cases of the solutions for the 
Whittaker-Ince limit of the GSWE and, for this reason, will be discussed in Sec. IV B. Here we 
regard the solutions W{u) for the Mathieu equation f|T2l) as a particular case (instead of a limiting 
case) of the solutions for the GSWE. Such solutions W{u) may be obtained by 

W{u) = U{z), z = cos^(^), (cr = l,z), 1 Mathieu equation 

} (48) 
zo = 1, 5i = -i, B2 = 1, 53 = 2P -a, uj = ±Ak, r] = } as GSWE, 

where U{z) are the solutions for GSWE given in Sec. II A, or, equivalently, by writing 

K = a/2, p= -1, ?9 = 4a - 8fc^ ^ = ±8ik (Mathieu equation as WHE) (49) 
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into the solutions W{u) for the WHE, given in Sec. II B. These solutions are also even or odd but 
their period are twice the period of the solutions corresponding to the WHE (if a = 1), that is, 
one has period 27r and 4n. As usual, there is no finite-series solutions. 
First set. Even solutions with period 27r if a = 1. 



cos^^ ) , |cosf^| > 1, (50a) 



W?{u) = e±4..cos^ (W2) g ^^i^j^^^cW sin^" |cos f | < oo 

oo . 
n=0 ^ ^ 

Wi{u) = e±4ifccos2(a«/2)^c(,^)|, ln + -,n + l]T8ikcos'^^j , Vu 

n=0 ^ ^ 

where the recurrence relations for the coefficients are 

{n + l)c% +[n{n± 8ik) ± 2ik + 2k'^ - a] c^^^ ± m - c^^l^ = (50b) 

In these and in the following solutions the upper or the lower sign must be taken throughout. 
Second set. Even solutions, period 47r if a = 1 

W^{u) = cos(a«/2) e±^^'=™^^('^"/2) J2 r[J+ (1/2)] ""^ T' 1'°' ^1 < 

W^2°°(m) = cos((tm/2) e±4ifccos2(W2) ^^(2) J^^sifccos^ ^) " , |cosf I > 1, (51a) 

n=0 

oo 

W'2(m) = cos((tm/2) e±^^'=^°'''('""/2) ^ c^f^* (^n + 1, n + 2; T^ik cos^ ^) , Vm 



n=0 

with the recurrence relations 

1 



{n + l)cf + 



n(n + l± 8ik) + - ± 2ik + 2k'^ - a 
^ 4 



clf) ± 8iA;n - clfl^ = 0. (51b) 



Note that in W^{u) we have eliminated the function \1/ (n + 1, n + 2; =i=8iA; cos^((T'u/2)) by using 
the integral representation [36] 

/"OO 

r(a)^(a,6;x) = / e-^*t''-^(l + t)''^"~icit, (3^6 > 3?a > 0) (52) 
along with the definite integral [3j 



oo 



re--cit = (53) 
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W^(u) = sm(au) e^^ifccos^wa) ^ ^' .^.. c^^ sm^" — , |cos ^| < oo 

' ^ \\n-\-i 



n=0 

oo 



Third set. Odd solutions, period 27r if cr = 1. 

r[n + (3/2)]^" 2 ' '^""2 

W3°°(m) = sin((TM) e±4^'^^°^'('^"/2) ^c^f)^ U + -,n + 3;T8zA;cos^ y) ' |cosf | > 1, (54a) 

Wz{u) = sin (aw) e±4^'^=°'^'('^"/2) ^ c(f)<| ( ^ + -, n + 3; T8ik cos^ — J , Vm 

n=0 ^ ^ 

having the recurrence relations 

(n + l)cjfji + [n (n + 2 ± 8ik) + l±6ik + 2k'^ - a] c|f) ± 8ik (n + cjf2i = 0. (54b) 
Fourth set. Odd solutions, period 47r if cr = 1. 

iy4°(«) = sin(aM/2) e±4^^^-'(-"/2) Y\n~^(3/2)f ^"^ ^' ^ 

„_n ~^ V / /J 



2 I ^ 



n=0 

-n-1 



H^4°°(m) = sin(aM/2) e±4^fc™^'('^«/2) ^ c^f^ (^^SiA; cos^ ^) " , |cos ^| > 1, (55a) 

n=0 

W^iu) = sm{au/2) e±4ifccos2{a„/2) ^ ^(4)$ ln+-^n + 2; ^Sikcos'^ — J , Vm 



n=0 

with 

1 



{n + l)citi + 



n (n + 1 ± 8ik) + - ± 6ik + 2k'^ - a 
4 



...,S*.(,..i)c-.^0. (55b, 



Finally, further solutions are obtained by performing the substitution 

71" 

u ^ u -\ — (56) 
a 

into the above solutions. This substitution, which leaves the Mathieu equation (1121) invariant, 
plays the role of the transformation T3T4, for it transforms z = cos^((Tm/2) into z = sin^(cr'u/2) 
without changing u. 

III. Solutions for the double-confluent Heun equation 

In the present section, first we obtain two sets of solutions for the DCHE 



z' 



+ {Bi + B2Z) -1- + (S3 - 2riujz + u'^z'^) U = 0, {Biy^O, uj ^ 0) 



dz'^ dz 



by taking the Leaver limit 2:0 of the first and third sets solutions for the GSWE given in Sec II 
A. After examining the general case and discussing finite-series solutions for a Schrodinger equa- 
tions with an asymmetric double-Morse potential, we establish the properties of these solutions 
for the Whittaker-Hill and Mathieu equations. 
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The following and f/g ^^6, in fact, the normal Thome solutions in the neighborhood of 2; = 0. 
As mentioned before, they encompass solutions which were investigated by Decareau, Maroni and 
Robert [Hj and also by Schmidt and Wolf [17] for a DCHE with four parameters. However, each of 
these solutions, resulting of the Baber-Hasse expansions, are now associated with two expansions 
in series of confluent hypergeometric functions and these converge in the vicinity of 2; = 00. 

A. General case and solutions of Decarreau et al. 

When zo —>■ 0, the limits of the solutions and Ui given in Eq. fl33al) are obtained immedi- 
ately. The limit of Ui yields the expression 



n=0 



where the recurrence relations for bn^ are 



Bi{n + + [n{n + B2-l)+ iuBi + ^3] b^n ^ + 2iu (^n + ir] + ^ - 1^ b';^!^ = 0. 



Then, by writing bn^ = [Bi) "'Cn \ we find the solutions given in Eqs. fl57ap and (I57b|) . 
First set. 

f/°(.) = e-J]cW(^-|-j , 

U^{z) = e^^^ ^ ci^)^ f n + zr/ + ^, n + B2] -2iujz j , (57a) 

n=0 ^ ^ 

^ f B \ 

U^{z) = e^^^^cW<l in + i7] + ^,n + B2,-2iujz] , 

n=0 ^ ^ 

where the recurrence for relations Cn ■* are {ci^\ = 0) 

(n + l)ci^i + [n{n + B2~l)+ luB^ + ^3] c« + 2iujBi (^ + ir, + ?^-l^ c^^\ = 0. (57b) 
Second set. To obtain the Leaver limit of the solutions given in Eqs. fl37al) we note that 



Hence, replacing 6^f^ by {—Bi)~'^Cn \ solutions fl37al) lead to 

00 / \ n 

U^{z) = e^^^+(-BiA)22-B2 ^ / ^ + 2 + ir] ^, n + 4 - S2, -2iujz , 

n=0 ^ ^ 

°° / B \ 

U2{z) = e*^^+(SiA)z2-fi2 ^ U + 2 + ii] -,n + A-B2, -2iujz 1 , 

n— n \ / 



(58a) 



with the recurrence relations 



(n + l)c% + (n + 3 - ^2) + 2 - luoB^ - + ^3] c^n^ 



~2iujBi in + 1 + ir] — 



B. 



0. 



(58b) 



From the previous solutions we may obtain others by the changes [rj^uj)^ {—Vy —^^)- 
Schrddinger equation with asymmetric double-Morse potential. Now we consider the time- 
independent Schrodinger equation for particle with mass m and energy E, 



du^ 



+ [£- V{u)]il) = 0, u:= Xx, £ 



2mE 



(59) 



where A is a real constant, x is the spatial coordinate and V{u) is a function proportional to the 
potential. For the asymmetric double-Morse potential considered by Zaslavskii and Ulyanov [HI [7] 
V{u) is given by 



/ (J 

V{u) = — I sinhw — — 



4 



B 



B \ s + -\ cosh u 



(60) 



where i? > 0, C > and s is a non-negative integer or half-integer. This is said to be a quasi- 
exactly solvable potential because a part of the energy spectrum and corresponding eigenfunctions 
can be found exactly in closed form |8l [9] . Here this means that one part of the energy spectrum 
is determined from the recurrence relations of an eigenfunction given by finite series. In section 
VI B, we show that the other portion can be obtained from the recurrence relations associated 
with an infinite-series eigenfunction. 

For the potential fl60l) the Schrodinger equation assumes the form 



+ 



4 8 16 



B fC 1 
-'2-2-^'^ 



B C 1 

— — H h s 

' 2 2 



du^ 

The substitutions [H [38] 

z = e", ij{u) = 0(z) = e-'''^^''^z^^'^^-'U{z), z e [0, 00) 
transform the previous equation into the DCHE (j4]) 



52 
16' 



^2u 



ij = o. 



(61) 



1^ 



C-2s)z 



dU 
dz 



' . 5' 2 



r. B (C \ , 
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t/ = 0, 



and then we can choose the following set of parameters 



5i 



B 



B, 



1 + C -2s, B-.=S 



B^ 
8 



-|- s — sC, iuj = — 



B 



IT] 



' 2" 



— s. 



(62) 



Hence, ^/'(m) = (l){z) is constructed by inserting into Eq. fl6Tl) solutions of the DCHE with the 
parameters specified in Eq. fl62|) and demanding that these eigenfuctions satisfy the regularity 
conditions 



lim ( 

2-^0 



0, lim 



0. 



(63) 
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Then, from the first set of solutions for the DCHE, we find 

oo 

n=Q 

oo 

n=0 

whose recurrence relations for bn^ are 

(n + l)b^^l^ +[n{n + C-2s)+£ + s^- sC] 6« - ^ {n - 2s ~ 1) h^^\ = 0, feL'J = 0. (65) 

Since s is a non-negative integer or half-integer, these are finite-series solutions for the reason given 
at the beginning of Section II. In effect, the coefficient 'jn^ of b^li vanishes for n = 2s + 1 and 
consequently the summation is restricted to the range < n < 2s. The regularity conditions are 
satisfied by both solutions because 5 > in the exponential factor; however only one solution is 
necessary to solve the problem (the series of regular hypergeometric functions is also appropriate). 
For a given value of s, the energies may be found by computing the determinant of the tridiagonal 
matrix (JTSj) corresponding to Eq. (!65l) . with < n < 2s. 

First, notice that the above solutions impose no additional restriction on the parameter C, 
contrary to the solutions employed in earlier approach pH]. Second, the solutions obtained from 
the other set of solutions for the DCHE are given by infinite-series but do not satisfy the regularity 
conditions: a suitable two-sided solution {—oo < n < oo) is presented in section VI. B. Third, for 
C = the potential becomes a symmetric double-Morse potential while the Schrodinger equation 
reduces to the Whittaker-Hill equation. In this case, by considering the WHE as a GSWE, we 
find regular infinite-series solutions given by one-sided series, that is, series with n >0 [TU] . 



f2z 



n 



2s, n 



C-2s- 



Bz 



(64) 



B. Whittaker-Hill equation as a DCHE 



Here we find that the solutions for the Whittaker-Hill, considered as a DCHE, possess no 
definite parity, contrary to the solutions which result from the GSWE (Sec. II B). Besides this, 
the periodicity (if any) will depend on the value of the parameter p. 

The substitutions 



z = e 



2iKu 



W{U) = ^l+(P/2)g5/(8.)^^^) 



transform the WHE 



dv? 



^ --i'^ - {p+ l)i cos(2/tM) + cos(4km) 



= 0, 



into 



,d'^U 
~d^ 



+ 



-| + (p + 3)z 



dU 
dz 



+ 



P 



+ +32-4 + 



-z z 
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U = 0. 



Then, by comparing this with the DCHE (j4]) we conclude that the solutions W{u) for the WHE 
are given by 



5i = -|, B2=p + 3, (| + l)Vg-f, ^.; = ±|, ^r^ = ±l{p+l) J DCHE, 
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(66) 



where now U (z) denotes solutions for the DCHE with the parameters specified in the second fine 
of the above expressions. 

First set. Taking ioj = ^/S and irj = {p + 1) /2 and using the solutions fl57ap . we find 



p^O(^y-) ^ g(g/4) cos{2Ku)+i{2+p)KU 



n=0 
oo 



n=0 ^ ^ 

oo . . 

n=0 ^ 



^2^nu^ > 0, (67a) 



where the recurrence relations Cn ^ are {cl\ = 0) 

(n+l)ci'ji+ n(n + p + 2)+ (l + l)'-^ 



.(1) _ i!. 



(67b) 



If p is a negative integer equal or less than —2, the series in these solutions are finite since the 
coefficient of c^-i vanishes for some value of n. For k = 1 the period of the solutions depends on 
the factor exp {i{p + 2)m). For instance, the period is: vr if p is an even integer; 27r if p is odd; 
2m7r (m > 1) if p = l/m, where / and m are prime to one another with / < m. 

Second set. Taking again iu = ^/8 and irj = (p + l)/2 but using the solutions (ISSal) . we 
find 



oo . s n 



n=0 

oo / 

W^{u) = e«/4)«i'i(2'^«)-^P'^« ^ c^f)^ ( n + 1, n + 1 

n=0 ^ 

oo ^ 

W2{u) = e«/4)Bm(2««)-ipK« ^ c(f)<l ( n + 1, n + 1 - 



7T _^p2iftu 



,2iK«| ^ Q 



(68a) 



p; 



n=0 
.(2) 



where the recurrence relations for Cn are 
+ l)cSi 



C(2) + i!!^c(2) _ 

16 ~ 



(68b) 



For these solutions the series are infinite. 

Third set. Taking iu = —^/8 and irj = — (p+ 1)/2 and using the solutions fl57ap . we find 



oo . . n 

W^{u) = e~(^?/4)sin(2Ku)+i(P+2)K« ^ / __ I ^(3)g2mK«^ 

n=0 ^ '^^ 



lyoo^^^ _ g-(i5/4)sin(2KM)+i(p+2)KM 



f2c^^^(n + l,n + 3+p;y'A 



= e-(^«/4)sin(2K«)-ipK« ^ c(f)<l> ( n + 1, n + 3 + p; ^e^*'"" 



n=0 



g2«K«| > 0^ ^gg^^ 
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where the recurrence relations for Cn^ are 



n(n + 2+p) + f-+ l + 



P 



+ -0 
16 ~ 



(69b) 



This set can also be obtained from the second one by p — >• —p — 2, ^ — > — ^. For these solutions 
the series are infinite. 

Fourth set. Taking iui = — ^/8 and ii] = — (p+l)/2 and using the solutions (I58ap . we get 



oo . s n 
-{$,/ A) cos{2 Ku)-ipKU \ ^ I _ \ J4) 2iKnu 

n=0 



?i=0 ^ 



-,2iKU I 



>0, 



(70a) 



Wi{u) = e-(«/^) cos{2^n)-ipK« ^ c(f)<l ( n - p, n + 1 - p; ^e^^'^'' 



n=0 



where the recurrence relations ch^ are 



n (n — p) + 



p V 
T ~ 4 



(70b) 



This set also follows from the first one by the changes C, and p — > —p — 2. If p is a positive 

integer, the series in these solutions are finite. 

C. Mathieu equation as a DCHE 

The solutions W{u) for the Mathieu equations (JT2i) are obtained from the solutions given in 
Sec. Ill B for the WHE by means of 



p 



■1, 2k = a, v = Aa- 8k^, ^ = ±8ik. 



(71) 



Unlike the solutions given in Sec. II. 3, now there is no solution with definite parity and for 
a = 1 {u real) the period is An instead of 27r and Air. We get rid of the functions which appear 
in and 14^2°^ by writing 



n+l 



^ (n + 1, n + 2; ±2^A;e'"") = (^±^ j £-^-("+1)-^ 

where we have used equations (152|) and (l53ll . 

First set. This comes from the first or fourth set of solutions for the WHE. 



(1) ianu 



n=0 ^ '^^ ' 

oo ^ 

y/^{u) = g±2ifccos((7«)-(i/2)(TM I J_ j (XK 

n=0 ^ ^ 

oo 

Wliu) = e±2ifc™«('^«)+(V2)<x« ^ c^i)!. (n + 1, n + 2; ^2ifce^''") , Vm, 



> 0, 



(72a) 



n=0 
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where the recurrence relations for ci}^ are (c*i{ = 0) 



{n + l)c^^l + 



n(n + l) -\ h 2A;^ - a 



Second set. This results from the second or third set of solutions for the WHE. 



W^{u) = e 



W^(u) 



q:2fesin(o-«) + (i/2)o-M ^ ^ / -p 
n=0 

CXD 

3=p2fc sin((Tu)— {i/2)<TU 

n=0 
oo 



2A; 



(2) io-nit 



1^2 (m) = eT2fcsin(an)+(i/2)a« ^ ^ $ (n + 1 , n + 2; T2iA;e^"") , Vm, 



< oo, 



> 0, 



n=0 

(2) 



where the recurrence relations for Cn' are [clI = 0) 



72 (n + 1) H 2A;^ - a 

4 



(72b) 



(73a) 



(73b) 



IV. Whittaker-Ince's limits for the GSWE 

Four sets of solutions for the Whittaker-Ince limit (IT!) of the GSWE, that is, for 

z{z - ^o)^ + {Bi + B2z)^ + [Bs + q{z - Zo)] U = 0, {q ^ 0), (74) 

result from the solutions given in Sec. II A by means of the very Whittaker-Ince limit {u — 0, 
7] oo such that 2riuj = —2q). However, we consider only the limit of the first set and, then, use 
transformation rules to obtain the other sets of solutions. In effect, if U{z) = U (Si, B2, B-^; Zq, q; z) 
denotes one solution for Eq. flT^ . the others are obtained by using the rules 

T^U{z) = z'+^^/^'>U{Cu C2, C3; Zo, q; z), Zo ^ 0, 

T2U{z) = {z~ z,f-^--B^I^^^U{B^, D2, D3; Zo, q; z), z^^ 0, (75) 

TiU{z) = U{-Bi - B2Zq,B2,B2, - qzo;zo, -q; zq - z), 

where Cj and Di are defined as section II. We again use only Ti and T2 but notice that T4 may be 
used to get alternative representations for the solutions. As a particular case we set up solutions 
for the Mathieu equation. 



A. General case 

First we rewrite the solutions given in Eq. (I33al) as 



00 



n=0 

00 

n=0 

00 



1 M ■ 



T I ■ B2 qz 

W \ n + iri + —,n + B2] — — , , 
2 ir]^ 



CXD X 

Ui (z) = e'^' ^ c^^^^l* ( n + ir] + 

n=0 ^ 



B2 , R qz 
2 irj 
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where we have multiphed the solution U^{z) by T ^1 + ii] — -y-] and written q = —2r]uj. From 
these solutions and respective recurrence relations for their coefficients, the limits are found by 
considering n fixed (a ~ it]) and using the formulae [35] 

lim $ fa, 6; = r{b)x^^-'^/^ Jh-ii2^), 

r(a + 1 - b)^ (^a, ^)] = 2x^^-^^/^Kb-ii2^) 



lim 

a— >oo 



where Ja(0 -^a(0 denote the Bessel and the modified Bessel function of the second kind 
whose definitions in terms of confluent hypergeometric functions are 

■Mt) = f(ATT) e-'* (I)' $ (A + 1, 2A + 1; 2^t) , ^^^^ 
Kx{t) = K_xit) = e-\2ty^ (A + i, 2A + 1, 2t) . 
For using (176|) when irj ^ oo {n and -B2 fixed, g =constant), ffist we make the approximations 

"if (n + iT] + ^,n + B2; -—) ^ ^ ( ir], n + B2; - — 
\ 2 ^^7 / V "^V 

^ ( n + if] + n + B2] -—] ^ $ fir], n + B2]- — 
\ 2 "^V J \ "^V 

In this manner we find the ffist set of solutions written in Eq. fl78ap . Notice that, although this 
is a formal derivation, the solutions can be tested by inserting them into Eq. flTll) . as shown in 
Appendix A. 

First set. These solutions admit the Leaver limit. If B2 + Bi/zq = 1, the ffist and the fourth 
sets of solutions are equal to one another. 

OD 
71=0 

00 

U^iz) = ^ c«(^y^)-"ir„+s,_i(2zy^), (78a) 

n=0 

00 

n=0 

where the recurrence relations for bn'' and c™ ■* are 

zo (^n + B2 + (n + l)b% + [n (n + S2 - 1) + ^3] 6« + qb^^l, = 0, (78b) 

and 

(n + l)ci'|i + [n (n + B2 - 1) + ^3] c^n^ + qzo (^n + B2 + ^ - 1^ c^^l, = 0. (78c) 

The coefficients bn^ and Cn^ are again connected by Eq. From this ffist set of solutions we 

obtain the others by applying the previous rules Ti and T2 according to relations (!32|) . 
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Second set. These solutions do not admit the Leaver hmit. If B2 + Bi/zq = 1, the second and 
the third sets are equal to each other. 

00 

n=0 

00 

U^{z) = ^(^-^^)/^5^c(fn^v^)-"^„+i+B.+H£i(2^v^), (79a) 

00 



n=0 

where the recurrence relations are given by 

Zo(n + B2 + ^^ {n + l)b% + 

2Bi\ f Bi 
n{n + l + B2 + — + + — 



^2 + — ) + ^3 



&f + = (79b) 



and 



[n + l)c2ii + 



n{ n + l + B2 + 



2Bi 

Zq 



^0 / \ Zq 



qzo \^n + B2 + ^-l^ c^^l, = 0. 

The relation between bn and Cn is given by Eq. ( !36l) . 
Third set. These solutions admit Leaver limit. 

00 

U',{z) = z'^^iz - zo)'-'''-^Y.^'i\z - zoT, 

n=0 

00 

U^{z) = z^^^iz - zo)'-'''-^Y.''n\^v^r''K^+3-B,{2tVQ^), 

n=0 

00 

U,{z) = z^^^{z - zo)'-""^^^ $^c(f)(^)-V„+3^B,(2v/^), 



n=0 



with the recurrence relations 

Zo 

and 

[n 



(79c) 



(80a) 



(^ + 2- B2-^{n + l)h% + [n{n + ?>- B2) + 2- B2 + B-^] hf + qh^^l^ = (80b) 



+ l)c% + [n{n + 3-B2) + 2-B2 + B^] c^f) + qzo (n + 1 - B2 - ^ j c^l, = 0. (80c) 
Eq. (138|1 affords the relations between bh. and ch . 
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Fourth set. There is no limit for zq ^ 0. 



UHz) 



UAz) 



,l-B2- 



[Z - Zo 



^ zo 2 2 



n=0 



[z - Zo 



oo 



£l+£2_l 



n=0 

OD 



n=0 



with the recurrence relations 
Zo f n 



Zo J 

Zo Zo \ 



Zo 



— -1+B, 



and 



in + l)citi + 



n[n + l — B2 — 



2Bi 



Zo 



Zo 



52 + — - 1 ) + 53 

Zo 



qzo[n + l-B2-^) = 0. 



(81a) 



(81b) 



(81c) 



Now the coefficients bn^ and Cn^ are connected by Eq. fHOj) . 

For the aforementioned wave equation which arises from the separation of variables of the 
Laplace-Beltrami operator for a scalar field on the Eguchi-Hanson space, Malmendier [19] has 
also found a solution in power series whose coefficients satisfy three-term recurrence relations. 
However, as a solution convergent in the neighborhood of 2; = 00, he has proposed a subnormal 
Thome solution with four-term recurrence relations, in contrast with the above solutions in series 
of Bessel functions. 



B. Mathieu equation and Lindemann-Stieltjes' solutions 

By inserting 

^0 = 1^ Bi = -1/2, B2 = l, z = cos^ (au), W{u) = U{z), 
into Eq. fm|) we obtain the Mathieu equation 



dv? 



+ [2q - AB3 - 2gcos(2(TM)] W = 0. 



Then, solutions for the Mathieu equation are obtained by making 

W{u) = U{z), z = cos^((T-u), (a = l,i), Mathieu Eq. as Whittaker- 

Zq = 1, 5i = -1/2, B2 = l, S3 = i[2g - a], q = A;^, J Ince limit of the GSWE, 



(82) 



where a = 1 for the Mathieu equation, a = i for the modified Mathieu equation and U{z) are the 
solutions for the Whittaker-Ince limit of the GSWE, Sec. IV A. Amongst the following solutions. 
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included in the W^{u), are solutions found by Lindmann and Stieltjes, which are even or odd and 
have period tt or 27r unlike the solutions of Sections II C and III C. As mentioned at the beginning 
of Sec II A, the following solutions may also be generated by letting ^ and p to tend respectively 
to zero and infinity so that pC, = 2fc^ in the solutions of Sec. II B to WHE ffTTj) {k = a), that is, 
by following the original proposal of Whittaker and Ince [21 |3l [23] . 

First set. Even solutions. The solutions Wi and Wi have period vr if a = 1 since Jn{—x) = 
(-irJJx). 



-i)"c; 



nJl) 



sin''" (ctm), 



cos (JM < OO, 



t^n-+(i/2)] 

OO 

W^{u) = ci^) [ik cos(orM)]"" Kn [2ik cos(aM)] , | cos(aM)| > 1, 



n=0 

OO 



Wi{u) = [fccos(aM)]"" Jn [2kcos{au)] , 



n=0 



The recurrence relations are 



(n + l)cl^l + + 



2 4 



+ k' 



n 



c^') -0 



Second set. Even solutions. If a = 1 W2 and W2 have period 27r. 



W^iu) = cos((tm) J2 



(1) 



?i=0 



r[n + (l/2)] 



sin'^fcTM) 



cos{(Ju)\ < 00, 



1^2°° (m) = ^c|,^^ [zfc cos((Tm)] [2iA;cos((TM)] , | cos(a?i)| > 1, 

n=0 

00 

W2{u) = ^c^^-* [/ccos(crn)]"" Jn+1 [2A; cos(crM)] , Vu, 



n=0 

with the recurrence relations 

.(2) 



{n + l)Ci + 



1 k"^ a 
n{n + l) + - + --- 



(2) 
n-1 



0. 



Third set. Odd solutions. If a = 1, and Wj, have period tt. 

00 { ^\n (3) 
H^30(«) = cos(aM) y -'.^.^ sin^^+ifcTM^ 



cos Cr-U < CX3, 



r[n+(l/2)] 

ly^(u) = tan(crM) Y^^n ^kco^ipu^"^ Kn^2 [2zA; cos(o"n)] , | cos(crn)| > 1, 

n=0 

00 

W'iiu) = tan((TM) c^^-* [A; cos((tm)]^" J„+2 [2/i; cos(crM)] , \\/u. 



n=0 

00 



n=0 



(83a) 



(83b) 



(84a) 



(84b) 



(85a) 
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The recurrence relations are 



(n + l)c% + 



/ „ \ k CI 

n n + 2 + 1 + — - - 
^ ^ 2 4 



.(3) 



Fourth set. Odd solutions. For a = 1, (z) and W^^z) have period 27r. 



1^4° M 



n=0 



r[n + (3/2)] 



sm 



2n+l, 



au), 



cos cr-u < oo, 



W^{u) = tan((TM) ^^c^"^"* [iA; cos(crn)] "iCn+i [2zA; cos(cru)] , | cos((jn)| > 1, 



n=0 

oo 



W4(m) = tan((TM) c^''-' [A; cos((Tm)] " J„+i [2A; cos(ctm) 



Vm. 



n=0 



The recurrence relations are 
(n + l)cSi + 



1 fc^ a' 
n(n+l) + - + --- 



From the preceding solutions other expansions may be obtained by the substitutions 

u ^ u -\ and k — > ik, 

2a 



(85b) 



(86a) 



(86b) 



(87) 



which plays the role of the transformation T4 given in Eq. 0751) . The solutions W^, W2 and 
those resulting from these by means of flHTl) . if expressed in terms of z [z = cos^ (c^m)], are the 
ones proposed by Lindemann and Stieltjes [291 EH] but here to each of these now correspond two 
solutions in series of Bessel functions. 



V. Whittaker-Ince's limits for the DCHE 

The solutions for the Whittaker-Ince limit of the DCHE, that is, for 

+ (5i + B^z)^ + {B, + qz)U = 0, {q ^ 0, B, ^ 0) (88) 

may be obtained from the solutions for the DCHE, given Sec. HI A, by letting u; — and t] ^ 00 
such that 2'quj = —q. However, at this stage, it is simpler to take the Leaver limit of the 

first and third sets of solutions given in Sec. IV A. As in the cases of the Mathieu equation there 
is no finite-series solution. 

First set. If we put bn^ = Cn V(-Si)" the limits of solutions fl78al) . we find 

00 / \ n 

00 

U^iz) = z(i-^^)/2^c«(^y^)-"ir„+5,_i(2zy^), (89a) 

n=0 

00 

Ui{z) = ^(i-^^)/2 5^c«(v/^)-"J„+s,_i(2y^), 

n=0 
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where the recurrence relations for c„ are (c_{ = 0) 

(n + l)cl'ji + [n{n + B2-l) + B^] c^^^ + g5ici'2i = 0. (89b) 
Second set. In the hmit of third set of solutions given in Sec. IV A we take 

Then, we obtain 

oo 
n=0 

oo 

U^{z) = ^(i-S2)/2gB,/. ^ c^^){^^)-K^_^,_^,{2^^), (90a) 

n=0 

oo 

f/2(z) = ^(i-^^)/2e^^/^^cf (^)-V„+3-i.,(2V^), 

n=0 

with the recurrence relations 

(n + l)c2ii + [n (n + 3 - 52) + 2 - ^2 + ^3] c^f) - g^ici'li = 0. (90b) 
As we see, in the present context, it is almost trivial to obtain solutions for Eq. (188|) . 

VI. Solutions with a phase parameter 

In this section we make some comments on solutions in which a phase parameter u was intro- 
duced to assure the series convergence if there is no free constant in the GSWE. Such solutions 
were partly studied elsewhere [H, [TOl [25l [26] but here we provide some more informations useful 
in applications. These solutions, which have been used in the study of the Teukolsky equations 
in Kerr spacetimes [25l EH ETJ [28], will as well be taken as a possible starting-point to find new 
solutions to the general Heun equation (Appendix B). 

Instead of three solutions as in Ref . |10j , we consider a basic set containing four solutions, two 
given by expansions in series of hypergeometric functions F{a, b; c; y) and two by expansions in 
series of confluent hypergeometric functions. These solutions, like the ones of Sec. II, admit both 
the Leaver and the Whittaker-Ince limits, but the last limit is not regarded here because it is the 
same as in Ref. [1], apart from the fact that now we have four solutions instead of two. 

In the first place we remark that there are several physical problems ruled by differential 
equations without any free constant, that is, requiring a phase parameter. For example, all the 
constants appearing in the GSWEs responsible by the time-dependence of Klein-Gordon and Dirac 
test-fields in some curved spacetimes are determined from the spatial part of the wavefunction 
[TOl [101 [m [12]. Similarly, in the radial Schrodinger equation for the scattering of electrons by a 
finite dipole [1] or by polarizable targets [U [ISj [HI [15] there is no disposable constant, since all 
the parameters which characterize the target and the incident particles are known. 

In the second place, solutions with a phase parameter u are given by doubly infinite series 
(—00 < n < 00) and, if 6„ denotes the series coefficients, then the recurrence relations have the 
general form 

anK+i + PnK + InK-i = 0, (-00 < u < oo) (91a) 
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where a„, 7„ and 6„ depend on v and on the parameters of the GSWE. These recurrence 
relations give the characteristic equation [1] 

a_i7o a_27-i a-37-2 , "oTi Qii72 «273 /^iun 

"» = ^aip P"') 

wherefrom must be determined. For a specific set of solutions we add a superscript in each of 
these quantities. 

Moreover, as we are going to explain, solutions with a phase parameter are also important 
because they can be used in two different manners to obtain solutions for equations having a free 
constant. 



A. Spheroidal equation 

The sets of solutions are denoted by by {U%, U°^, U^, U°^) [i = 1, 2): the solutions [7°, converge 
for any finite z and the others converge for \z\ > \zq\. These solutions were obtained by integrating 
the the GSWE [H [TOl [25l [26], except which is derived from solutions of the general Heun 
equation in Appendix B. We write only the first set of solutions (the second set is obtained from 
the first by using the transformation rule T2), namely, 

UUz)=e- b^^^F(^-n-u~l,n + u + ^;B, + ^;'-^), 

\ / Z Zq Zq J 

n=—oo ^ <j <j / 

Br, 

n=— 00 ^ 

F(n + u + f,n+u + l-f-^;2n + 2u + 2; ^) , (92a) 

00 

U^{z) = e^'^^zi-{^2/2) J2 &i^H-22CJz)"+'^^(n + u + 1 + ir],2n + 2u + 2; -2iujz), 



-n—v — ^ 



n=— 00 
00 



U^{z) = e^'^^zi-(«2/2) 6W(2icuz)"+^<l>(n + u + l + ir],2n + 2u + 2] -2iujz), 

n=—oo 

where the functions F{a, b; c; y) and $(a, b; y) are defined by means of 

Fia,b;c;y) = ^^^^, 6; ^) = £(iz^$(a, 6; ^) (92b) 

r(c) r(6) 

and the coefficients bn^ are connected with bn^ through 

6« := (-l)"r (^n + u + 2-^y(^n + u + l-^-^yil\ if Zo 7^ 0. (92c) 

The coefficients of the recurrence relations (19 lap for bn^ are 

(I) (n+u+2~-^)(n+u+l-^-^)(n+u+l-iv) 
an = iuJZn 7 -7 TT , 

=-B.- Vu:z, + . + 1 - f + . + f) - ^^tt:lt;^) , (92d) 



(1) . {n+u+^-l){n+u+?^ + ^){n+v+i^) 

In — —IOJZq -7— IN . . 

2(^n+!/-ij(n+!/) 
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Note that the hypergeometric functions in U^^^ are well defined only if B2 + {Bi/zo) 7^ 0, — 1, — 2, • • • , 
while in f/gj, we would have the restriction B2 + (Bi/zq) 7^ 2, 3, 4, ■ ■ ■ . 

In these solutions, u cannot be integer or half-integer in order to prevent two dependent terms 
(for different values of n) inside the summations. Suppose, on the contrary, that 2i/ + 1 = iV in 
for some integer A^. Then, for n = rii (Vni) and n = —N — rii, the hypergeometric functions 
in Ui^, would be equal on account of their symmetry with respect to the two first parameters. One 
may show that this is the same for the other solutions. In fact, integer or half-integer values for 
u would lead to infinity coefficients in the recurrence relations for bn\ 

Now we consider the two possibilities for getting solutions without phase parameter from the 
ones with a phase parameter. The first consists in breaking off the series by taking n > in which 
case the parameter u is given in terms of the parameters of the differential equation [U [TD]. In 
particular, the truncation of gives an expansion in series of Jacobi polynomials which includes 
the Baber-Hasse solutions (associated Legendre polynomials) for the angular equation of the two- 
center problem of quantum mechanics |24j as well as the Fackerell-Crossman solutions for angular 
Teukolsky equations [46j. In this regard, we must do the following correction in Ref. [lOj: the 
truncated expansions in series of regular confluent hypergeometric functions converge for \z\ > \zo\, 
rather than for any z as stated there. Therefore, such solutions are not suitable to solve an angular 
equation because in this case \z\ < \zq\. 

The other possibility gives solutions in terms of doubly inflnite series (—00 < n < 00). These 
are obtained by ascribing any convenient value to u, excepting integer or half-integer values. Such 
solutions exhibit properties that may be quite different from the properties of the one-sided series 
{n > 0), as illustrated by the following solution for WHE which we derive from a solution for the 
GSWE. Proceeding as in Sec. II B and using the relation [36] 

F(— a, a; 1/2; sin^ m) = cos(2au), (93) 

we flnd that the solution U^^^iz) affords an even solution for the WHE, namely, 

00 

Wliu) = e(«/2)-^'(-") J2 bn cos[(2n + 2u + 1)ku], (94) 

n=— 00 

where the recurrence relations for bn are (— 00 < n < 00) 

e (^n + z/ + ^ - bn+i + [^9 - (2n + 2u+ if] 6„ - ^ (^n + i/ + ^ + b^-i = 0. (95) 

Now, if there is a free parameter in the differential equation, we can choose any admissible value 
for u and, in particular, 

2z/ + l = //m, K=l, (96) 

where / and m are integers prime to one another, so that 1 < I < m. Then, we flnd the following 
two-sided Ince solution, W({u), for the WHE ^ 



(97) 



having period 2iTm, m > 1, in opposition to the solutions obtained from the truncated solutions 
[To] (period vr or 2n). An odd solution may also be found from the f/2iy(^) of the (omitted) second 
set of solutions. 
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Although this example is concerned only with the periodicity and parity of solutions for the 
Whittaker-Hill equation, the analysis can be extended to other aspects of the GSWE and their 
limiting cases. Next we use two-sided solutions similar to the above one to complete the solutions 
of the Schrodinger equation considered in section III. 



B. Double-confluent Heun and Schrodinger equations 

From the previous solutions there is no difficulty in finding solutions for the DCHE by means 
of the Leaver limit, provided that we use the formulas [35] 



c— >oo 



y 



lim F ( a, 6; c; 1 - - ) = y"^(a, a + 1 - 6; y). 



6— >oo 

We find the Leaver solutions [H fTO] 



B2 

\v 

n=—oo 



/ R \ "^'^ / R R \ 

E ^i'Mf) ^U + - + f,2n + 2. + 2;^), 

=— oo \ / \ / 



n=— OO 

OO 



V^^iz) = e'^'z'^-^ E h^^\-2iujzY^''^>{n + u+l + iT],2n + 2u + 2;-2iujz), 

n=— oo 
oo 

U'^{z) = e^'^^z^-^ E 6(,^)(2ituz)"+^l'(n + z/ + l+ir/,2n + 2i^ + 2;-2icuz), 



with the following coefficients in the recurrence relations fl91ap for hn^: 

(1) . (n+v+2-^)(n+u+l-i'q) 
an = lUjBi ^ ' 



(98) 



(99) 



4" = B,-\ (B, - If + 1)^ + 5;SlM. (100) 



(1) _ (n+u+^~l)in+u+iri) 

— LLUlJi ; -7 p> . 

' ^ 2(n+i/)(n+i/-|) 

These four solutions result as the Leaver limits of solutions for the GSWE because in Sec. VI. A 
we have chosen a set formed by four solutions. They can also be obtained by starting with the 
two expansions in series of confiuent hypergeometric functions for the GSWE and, then, using a 
symmetry of the DCHE [H [TO]. 

Now we use these solutions of the DCHE to find infinite-series solutions for the quasi-exactly 
solvable double-Morse potential discussed in section HI. For this is necessary to match two solutions 
having different regions of convergence, a procedure which has been employed to find solutions for 
the radial Teukolsky equations in Kerr backgrounds [25l [26] [27] [28] . In this manner, the missing 
portion of the energy spectrum may be obtained. 
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By inserting the previous 11^^^ and into Eq. (1611) and using the parameters given in Eq. 
f l62|) . we get, respectively, 



0U^)=exp[-f (. + !)] Y: ^i^Mi 

n=—oo ^ 

0SW=exp[-f (. + !)] f; 



1 C -B \ 

^(n + i/ + - + --s,2n + 2z/ + 2;— j. 



\ n+u+4, 

Bz ^ 2 



(101) 



1 c 



Bz 



^ {n + u + ----s,2n + 2u + 2- — , , 



where in the recurrence relations (19 lap for now we have {B > 0, C > 0) 

(l) ^ ^2 (n+^+|-g+3)(n+r/+| + g+^) 

nil) _C,Bi_C!,/^,^,n2_ [C^-(l+2s) 
Pn - C + g 4 + V'' + ^ + 2;' 32 (n+j.)(n+i.+ 



(n+j')(n+i/+l) ' 



(102) 



(1) ^ B2 (n+v+g-l-s)(n+z.-g-l-.) 



For infinite series, these solutions satisfy the regularity conditions fl63l) since we have {B > 0) 



lim 0ij^(-2) ~ 2; 2 exp 
2— >o 



B 



2; + 



2; 



0, 



lim 



6^(2) ~ 2 2+^ exp 



5 / 1 
4 V z 



0, 



where we have used the relation (1271) . However, in order to guarantee that in fact these are 
solutions given by infinite series, we must choose suitable values for the parameter u. Recalling 
that u cannot be integer or half-integer and that s is a non-negative integer or half-integer, we 
find that there are three different cases to be considered. 

First case: C is integer or half-integer. In this case we can take z/ = 1/3 to find 



16 



(n+l+i)(r 



3^2, 



(1) 



1: 



(1) 



P+B!_^^(„^l + 1^2_ [C^-{l+2sf] 
C+g 4+1^+3+2; 32 

■)(«- 



(103) 



16 



1 , 1 _ 



(„+l)(„+l_l) 



Thence, the numerators of ctn^ and '-)n^ do not vanish since (1/3) + (C/2) and (1/3) — (C/2) are 
not integer or half-integer. Therefore, we have a pair of solutions given by two-sided infinite series. 
Second case: C is not integer or half-integer, but s is integer. In this case we can take 

C , 

which is different of integer or half-integer. Then, 



In 



0("+^ 



16 („+2+f)(n+| + f) ' 

_ ('^-^+i+C)("-^+l) 
16 ^n+l+%)(n+l + %) ■ 



B2 [c2-(1+2s)^] 
32 („+i+C)(„+2+£) 



(104) 
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Thus the numerators of and "jn^ cannot vanish and, once more, we have two-sided infinite 
series. 

Third case: C is not integer or half-integer, but s is half-integer. In this case we can take 

C 1 

V = h - 

2 2 

which again imply that the numerators of ai^^ and 7!^'' cannot vanish since 

(1) _ b2 {n+s+2){n+s+2+C) 
- 16(„+3 + C)(„+2+f)' 

^(1) _ C I i?^ \ in \ l \ \C^-{l+2s?] .^Q,^ 



2 ' 2 , 



^,(1) _ {n-s+C){n-s) 



Therefore, in principle this problem can be solved by matching the two double-sided solutions 
which converge over different ranges of the independent variable z. The energy spectrum can be 
obtained by solving the characteristic equation fl91bp . Then, for a given value of s, besides the 
energy levels provided by the quasi-polynomials solutions of Sec. III.A, there is the spectrum 
resulting from the present infinite-series solutions. In addition, for C = 0, the first case affords 
solutions in terms of doubly infinite series for the symmetric double-Morse potential (WHE). 



VII. Conclusions 

We have seen that the Leaver and Whittaker-Ince limits of the generalized spheroidal wave 
equation (GSWE) allow us to generate solutions for three other equations with different types of 
singular points. We have also obtained new solutions for the Mathieu and Whittaker-Hill equations 
by regarding these as particular cases of the GSWE and its limiting cases. As an application we 
have considered the Schrodinger equation with an asymmetric double-Morse potential. 

First we have constructed a set of solutions to the GSWE in such a manner that the three 
solutions of each set converge over different regions but have series coefficients proportional to 
one another. Actually, this is an extention of an approach already used in a previous work [lOj 
(revisited in Sec. VI), the difference being that in Sec. II the initial set of solutions is given by 
a series of ascending powers of 2; — z^^ due to Barber and Hasse, and two expansions in series of 
confluent hypergeometric functions. From the initial set other ones have been generated by means 
of transformation rules which reflect the symmetries of the GSWE. The importance of these rules 
can be appreciated by considering the properties of each set of solutions that we have found for 
the Mathieu and Whittaker-Hill equations. 

From solutions of the GSWE, sets of solutions for the three limiting cases have been found 
following the diagram ffTOl) . Thus, the solutions for the double-confluent Heun equation (DCHE) 
have resulted by means of the Leaver limit. On the other hand, in the Whittaker-Ince limit for the 
GSWE and DCHE, the expansions in series of confluent hypergeometric functions have afforded 
expansions in series of Bessel functions as solutions for Eqs. ([7]) and ([9]). 

In this process for generating solutions, the Barber-Hasse expansions have led to some known 
solutions for the DCHE, Mathieu and Whittaker-Hill equations, but these solutions are accom- 
panied by two new solutions belonging to each set. Besides this, the Whittaker-Ince limit of the 
Barber-Hasse expansions has also given solutions which include as particular case the solution 
already proposed by Malmendier [19] to solve a wave equation which arises from separation of 
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variables of the Laplace-Beltrami operator for a scalar field on the Eguchi-Hanson space. How- 
ever, now we have found two expansions in series of Bessel functions, both of them convergent in 
the vicinity of 2; = oo and satisfying three-term recurrence relations for the series coefficients, in 
contrast with the solution proposed by Malmendier (four-term recurrence relations). 

The properties of the solutions for the Mathieu and Whittaker-Hill equations depend on 
whether these solutions come from the GSWE or DCHE. We have seen that the unique solu- 
tions for the WHE which are even or odd and have period vr or 2tt result from the solutions 
for the GSWE (Sec. II B). The Whittaker-Ince limit transfer these properties to the solutions 
of the Mathieu equation given in Sec. IV B, each set including one of the Lindemann-Stieltjes 
solutions and two new solutions in series of Bessel functions. The latter solutions demand further 
investigation since they are different from the usual expansions in series of Bessel functions. 

On the other hand, the solutions for WHE considered as a DCHE (Sec. HI B) do not possess 
such properties of parity and periodicity (even or odd, period vr or 27r). The same is true for the 
solutions of the Mathieu equation when this last is treated as a particular case (not as Whittaker- 
Ince limit) of the WHE (Sec. II C) or as a particular case of the DCHE (Sec. HI C). However, 
solutions for the Mathieu and Whittaker-Hill equations which do not present the usual properties 
of parity and periodicity may be important to solve initial value problems instead of boundary 
value problems [22]. 

Up to section V we have assumed that there is a free constant in the differential equation, and 
all the solutions are one-sided in the sense that the summation index n runs from to cxd. In 
this case, the free constant is determined from a characteristic equation which must be fulfilled 
in order to assure the series convergence. However, in section VI we have reconsidered a set of 
solutions in series of hypergeometric and confluent hypergeometric functions which converge when 
there is no free constant in the differential equation. In such solutions there is a parameter u, 
called phase or characteristic parameter [U |2T], which must be adjusted so that a characteristic 
equation is satisfied in order to guarantee the series convergence. 

The solutions with the phase parameter u are two-sided (—00 < n < 00) and also admit both 
the Leaver and Whittaker-Ince limits. In addition, they allow us to find one-sided and two-sided 
solutions for an equation having a free parameter. The one-sided solutions are obtained by the 
restriction n > in which case u becomes determined in terms of other parameters of the equation 
[TU] . The two-sided solutions are generated by choosing some of the admissible values for u, that 
is, some u different from integer or half-integer. These two-sided solutions have been used in 
section VI B to find infinite-series solutions for the Schrodinger equation with an asymmetric 
double-Morse potential, by reducing that equation to a DCHE. 

The asymmetric double-Morse potential considered here is a quasi-exactly solvable (QES) 
potential for which a part of the energy spectrum can be determined from the recurrence relations 
for the coefficients of finite-series eigensolutions, as the ones found in Sec II A. On the other hand, 
in Sec. VI B we have shown that, for a given value of the the parameter s, the other portion of the 
spectrum can be derived as solutions for a characteristic equation associated with eigenfunctions 
given by infinite two-sided series, provided that we accept to match two solutions which converge 
over different ranges of the independent variable. When the parameter C vanishes, the asymmetric 
potential degenerate to a symmetric one, whereas the Schrodinger equation degenerate to a WHE 
which has already been solved as a GSWE [TU]. Therefore, this symmetric case is an instance of 
a WHE which can be regarded both as a DCHE and a GSWE. 

On the other hand, for the trigonometric and hyperbolic QES potentials given by Usheveridze 
[8], the Schrodinger equation can be transformed into GSWEs [38]. Thus, the solutions for the 
GSWE given in this paper may be used to verify if it is possible to find the whole energy spectra 
for these potentials, as in the case of the double-Morse potentials. 
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Another possible application consists in taking the solutions for the GSWE as a guide to find 
new solutions for the (general) Heun equation, since the former equation is a confiuent case of 
the latter. This would complement the procedure presented in the preceding sections, because 
eventually it would add the connection: general Heun equation — > GSWE into the diagram ffTOl) . 
In fact, this link is required by the study of Teukolsky equations for gravitational backgrounds of 
black holes, as we have mentioned in the first section. Although the Heun equation is not subject 
of the present paper, in Appendix B we have introduced some preliminary results which indicate 
how the above goal could be attained. 
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Appendix A: Solutions in series of Bessel functions 

To find directly the solution U^{z) for the Whittaker-Ince limit of the generalized spheroidal 
wave equation. 



dz^' 



we perform the substitutions 



to obtain 



,d^Y dY 



{B, + B^z)^ + [53 + q{z - Zo)] U = 0, {q^ 0), 

t = 2i^, Ur,{z)=t'~''^Y{t), 

d^Y Aq{zQ - 2Bi - 2B2Z0) dY 



(Al) 



(A2) 



+ t— VY = -AqzQ- 



+ 



df^ dt ^ " df^ 

2Bi + ^22:0 + Zo 



4g(l-52)- 



+ (1 - B2Y + Aqzo - 4S3 



dt 
Y. 



Next, by expanding Y{t) in series of modified Bessel functions of the second kind K\{t) according 
to 



Y{t) = J2 dJ-'^Kxit), where X = n + B2 - I, 



(A3) 



n=0 



we find 

dY 
'dt 



n=0 



dK\ n 



dt t 



(PY_ 
' IF 



^dj 

n=0 



d'^Kx 2ndKx . n{n + l) 



dt^ 



t dt 



+ 



t' 



(A4) 



Inserting ( ]A3[) and (1A4P into the differential equation for Y(t) and using the modified Bessel 
equation. 



,d'^Kx dKx 



dt^ 



+ t- 



dt 
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it results 



J2 ^nt-" [A^ + n^ + 4^3 - Aqzo - (1 - ^2)'] Kx 



n=0 

00 



n=0 



rf^fsTA dKx , ,1 rf/^A 

-4g2o-5i^ + 2^^-17^ + 4g (2nzo - + 220^2 + 25i^ 



-4g^ ci„t-"-2 [n(zon + 2zo + 25i) - (1 - B2){zoB2 + 25i + z^)] K, 

n=0 

From the difference-differential relations for K\ l39i, we find that 



(fK, 

dt^ 
dKx 



(A5) 



t 



dt 



A 

' 7rA + A(A-l)t-2irA + t-'i^A+i, 

-Ai^A-t^A-l, 



Ar^irA-r^iTA+i. 



1^ 

t dt 



Thence, taking into account that X = n + B2 — 1, Eq. flASP simplifies to 



J2 5(^nr"+ifr„+B2-2 + J2^n{n + B2-1) + B3]dJ-^Kn+B,~i + 



n=l 
00 

E 

n=0 



n=0 



22:05' ( n + -B2 



Bi 



dJ-'^-'K. 



n+B2 



since the coefficient of t " ^Kx vanishes. By the replacements n — >■ m + 1 and n — m — 1 in the 
first and third terms, respectively, this equation becomes 

00 . ^ 

[Idi + Bsdo] Kb,-i + ^2)2^'' + i)dn+i+ [n{n + B2 - 1) + 53]rf„+ 

n=l 

2zoq (n + B2 + ^-l^ d„-i] r"i^„+B,_i = 0. 

Requiring that the coefficient of each independent term of this equation to vanish, we get the 
relations = 0) 

^(n + 1)4+1 + [n{n + B2 - I) + B^ld^, + 2zoq (n + B2 + — -l) d^-i =0, n > 0. (A6) 
2 \ zo J 

Then, from Eqs. flA2l) and flASp we have 

00 00 
UrA4t)) = t'~^^ 5Zf/„t-"i^n+B,-i,W oc 5^2""rf„(zy^)-"ir„+B,-i,(2^v^). (A7) 



n=0 



n=0 



By putting Cn^ = 2~^dn in Eqs. (1A6I) and ( 1A7I1 . we recover the solution given in Sec. IV. 
The derivation of the expansion in series of Jx{2^/qz) is similar to the previous one. Thus, the 
Whittaker-Ince limits found in Sec. IV in fact satisfy the differential equation flAip . 
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Appendix B: General and confluent Heun equations 



The relationship between solutions of general Heun equation and GSWE (confluent Heun) is 
interesting by itself and becomes, in the context of relativistic astrophysics, important in the study 
of the Teukolsky equations of the gravitational backgrounds of black holes. In effect, for the Kerr 
metric these equations turn out to be GSWEs which afford DCHEs in the extreme upper limit 
of the rotation parameter [Tj. On the other hand, the Teukolsky equations for the Kerr-de Sitter 
metric are (general) Heun equations which, in turn, reduce to GSWEs when the cosmological 
constant A vanishes [TTl [T2] . Therefore, in such problems we have the connection: Heun equation 

GSWE — s> DCHE, and this requires a generalization of the diagram ( ITOj) . 

In the following we exhibit two Erdelyi's solutions for the Heun equation [17] which lead to the 
expansions in hypergeometric functions given Sec. VI for the GSWE. From these Erdelyi solutions 
we explain how to obtain two additional solutions which yield the expansions in series of confluent 
hypergeometric functions for the GSWE as well. 

To begin with, let there be the Heun equation [151 HHl HS] 

+ - + 7 + + ^TTT-^ 7^ = 0, e = a + /9 + l-7-5, (Bl) 



dx'^ \x x — 1 X — a J dx x{x — l){x — a) 

where a G C \ {0, 1} and x = 0, 1, a, cxd are regular singular points with indicial exponents given 
by {0, 1 — 7}, {0, 1 — 6}, {0, 1 — e} and {a, /?}, respectively. For a = and a = 1 the equation can 
be reduced to the hypergeometric equation by changes of variables; also for e = and q = aa[3 
it degenerates to the hypergeometric equation. There is no possibility of confusing the auxiliary 
parameter q with the parameter q of the Whittaker-Ince limits given in Eq. ([6]) neither the 
singularity parameter a with the constant a of the Mathieu equation (fT2i) because the q and a of 
the Heun equation disappear in the confluent limit of the equation (see Eq. (1B6P ). 

As in Ref. [12], we denote by if|°'^^(a;) and h[1'^^ {x) the two Erdelyi solutions in series 
of hypergeometric functions, where the superscripts mean that the first solution converges in a 
domain containing the singular points and 1, and the second, in a domain containing a and 00. 
These solutions are [17] (once more v is different from integer or half-integer) 



00 / 



7+5 7+5 
-n-v -1 + , n + v + ; 7; x 



n=— 00 
00 



[B2) 



where 



6(1) = (-l)«r {n + iy + 2- F f n + z/ + 1 + ^-^ ) 6«. (B3) 
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In the recurrence relations (19 lap for bn^ the coefficients are 

(n + iy + 2-^^){n + iy+l-a + ^^){n + 1^ + 1-/3+ ^^){n + iy + l-^) 



4 (n + z/ + 1) (n + z/ + I) 



{l'"") ^""^"^ (^ + ^ + l) + ^[2«/?+(« + /?)(7-'5) + (7 + 5)(l-7)] 

1 (^_^)(^ + ^_2)(2«-7-^)(2/3-7-^) 
--a(7 + 5)(2-7-5)-g+ 32(n + .)(n + . + l) ' 

^(1) _ jn + u + a-^^) (n + u + P-^^) {n + u - 1 + {n + u + --^) 

4 + z/ - ij (n + z/) 

In these solutions the parameter u which appears in some other works, for example in Ref. 
|15j . have been replaced by the expression on right-hand side of the relation 

+ (B5) 

This redefinition is due to Suzuki, Takasugi and Umetsu [12]. It is useful to get the limits for the 
GSWE, as well as to obtain new solutions from the preceding ones by means of transformations 
like that written in Eq. (IBlip . simply because it is not necessary to transform u in order to get 
solutions having the same series coefficients. 

The above solutions are valid only if a ^ [0, 1] [47j. Since this restriction comes from the analysis 
of the recurrence relations for the coefficients bn^ as n — » ±oo, it is transferred to solutions whose 
coefficients are proportional to bn^ . The solution hI^'^\x) converges in the interior of an ellipse 
having foci at and 1 and passing through the point a, provided that the characteristic equation 
is satisfied; the solution hI'^'°°\x) converges in the entire complex plane excepting the line joining 
the points to 1 [47J. Therefore, we are dealing with a pair of solutions convergent over different 
domains. 

To derive the GSWE as a confluence of Heun's equation fIBip . ffist we rewrite the latter as 



X \ d'H 



Then, by letting that 



X A ^ (X \ t , dH 



7(s - 1) 1 + 5s 1 + -xix - 1) 

a J \a ) a 



dx ^ 



(3 q 

a—x 

a a 



H = 0. 



a, /?, q oo such that — — > > — p, - —a, H{x) h{x), (B6) 

a a a 

p and a being constants, we find the GSWE [15j 

d'^h dh 
x{x - 1)— + [-7 + (7 + 6)x + px{x - 1)]— + [apx - a]h = 0, (B7) 
dx'' dx 

which, however, is not in the form suitable to get the DCHE by means of the Leaver limit. Despite 
this, from the previous solutions, we see that H{x) and its limit h{x) have the same functional 
form. To obtain the recurrence relations for the series coefficients, ffist we divide the coefficients 
( 1B4|) by a and then take the limits of an^ /a, (3n^ /a and 'yn^ /a when a — 00, by using the relations 
(lB6p . To permit the Leaver limit, the second step is given by the substitutions 

X = h{x) = e''^/('^°)f/(z) (B8) 

Zq 



39 



which convert Eq. (lB7p into 



dU 



z{z -zo)-— + [Szo + (7 + S)z]— + 



dz 



dz 



IP 



Then, comparing this with the GSWE ([T]), we find that 



a = iri + —,d = , 7 = ^2 

2 zo 



_ 2 




H), 


we 


Bi 






P 







a — \ a 



7 + 5\ p{z- Zo) p^z{z - Zq) 



Zo 



^Zq 



U = 0. 



-2iujzo, a = —iuzQ ( B2 + 



Bi 

Zo 



(B9) 



Therefore, thanks to Eq. (IBSp . the solutions U{z) for the GSWE ([T]) are obtained by writing 

Zo - z 



U{Bi,B2,B3;zo,u,ri;z) = e"^^h a^-f, 6, p, a; x 



(BIO) 



and by inserting the parameters (IB9p into the right-hand side and also into the recurrence relations 
obtained in the first step. 

By this procedure, the solutions hI^'^\x) and hI'^'°°\x) for the Heun equation give, respec- 
tively, the solutions Ui^,{z) and U^{z) for the GSWE, written in Sec. VI . To obtain solutions 
which lead to the expansions U^{z) and U^{z) for the GSWE, we achieve an appropriate substitu- 
tion of the independent variable, followed by a transformation of the dependent variable, such that 
the Heun equation is transformed into another version of itself (with a different set of parameters). 
More precisely, from the transformation theory for the Heun equation ^9], if H{a, q; a, /3, 7, 6; x) 
is a solution of the Heun equation, then another solution is obtained by means of (transformation 
number 105 of Maier's table j49j ) 



H{a,q;a,P,j,S;x) 

{x — 1)~" H i a,q + a{a — 7 — 5 + l)a; a, a — 5 + 1, e, a — /3 + 1 



X — a 
X — 1 



(Bll) 



Then, by identifying if (a, q; a, (3, 7, 6; x) with hI'^'^\x) and if^"'°°^(x) we find that 



r{a,oo}. 



00 /■ 

'n——r^\ ^ 



7 + 5 . 7 + 5 X — a 

^,n + . + l + a-^;e;^ 



H[r\x) 



\x 



1)-" E 



1 - ^ 

a_ 

X — 1 



X 



(B12) 



F + z/ + 1 + a - n + z/ + 1 - /3 - 2±^; 2n + 2i/ + 2; fEi) 



where 



.(1) 



a)"r [n + v + 1- a- 



7 + 5 



V [n + u+l- (3 + 



Note that u has not been transformed and that the coefficients hn^ are the same which appear in 
the original solutions. 

Following the procedure outline above for a — *^ and using the limits fl98l) for the hypergeo- 
metric functions F(a, 6; c; y), we can show that these give the solutions U^{z) and U^{z) in series 



5 + 7 



(B13) 
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of confluent hypergeometric functions for the GSWE. In the first solution, before taking the limit 
of the hypergeometric functions we make the approximations 

X — a 1 — a c. 
c = e = —pa, = 1 H ~ 1 



X — 1 X — 1 p(l — x) 

In the second solution, first we must take 

X — 1 p(l — x) 



b ~ —P = pa, 

X — a 



and, in addition, we must find the recurrence relations for the coefficients Cn ■* defined by Eq. flB13l) . 
However, to establish precisely this first set of solutions, it is necessary to study the convergence 
of the two last solutions and, then, apply the so-called elementary power transformations of the 
dependent variable ^2] in order to generate new sets of solutions. We have also to find solutions 
valid for a e (0, 1). 

Finally we mention that the following series solution for the Heun equation [T5] 

oo 

= ^rf„x", |a|>l, 

n=0 

where the coefficients satisfy the relations 

a{n + l)(n + 7)c?„+i — [an{n + 7 + 5 — 1) + n{n + a + P — 6) + q]dn 

+ {n + a-l){n + p- = 0, 

converges for |x| < \a\ and gives, by the above procedure, the solution Ui{z) of section II for 
the GSWE. By analogy with the previous case, we may expect that the generalization of the 
solutions U^{z) and Ui{z) (in series of confiuent hypergeometric functions) are given by series of 
hypergeometric functions . 
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